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Abstract 

The well-known energy problem is discussed in f{R) theory of grav- 
ity. We use the generalized Landau-Lifshitz energy-momentum com- 
plex in the framework of metric f{R) gravity to evaluate the energy 
density of plane symmetric solutions for some general f{R) models. 
In particular, this quantity is found for some popular choices of f{R) 
models. The constant scalar curvature condition and the stability con- 
dition for these models are also discussed. Further, we investigate the 
energy distribution of cosmic string spacetime. 

Keywords: f{R) gravity. Generalized Landau-Lifshitz EMC, Energy 
Density. 

1 Introduction 

The energy locaUzation has been a thorny problem since the Einstein era. 
Several attempts have been made to find a general and unique tensor rep- 
resentation for the energy-momentum. Einstein was the first who tried to 
solve this problem by introducing energy-momentum pseudo tensors. He 
established the energy-momentum conservation laws given by [1] 

^-{v^(T; + g} = 0, (/i,i/ = 0,l,2,3). 



dx 
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where is the energy-momentum density of matter and represents the 
energy-momentum density of gravitation. It is mentioned here that tj^ is not 
a tensor quantity rather it is the gravitational field pseudo-tensor. Komar |2] 
gave a set of energy-momentum covariant conservation laws and developed 
their relationship to the generators of infinitesimal coordinate transforma- 
tions. Bergmann [5]- [5] contributed greatly to the fundamental nature of 
conservation laws. 

Landau-Lifshitz introduced [6] the energy-momentum complex by using 
the geodesic coordinate system at some particular point of space. Many other 
people like Tolman [7], Papapetrou [8], Bergmann [9], Goldberg [10], M0ller 
[TT] and Weinberg [12] developed their own energy-momentum complexes. 
All these prescriptions, except M0ller, are restricted to perform calculations 
in Cartesian coordinates only. Also, we cannot define angular momentum 
with the help of these prescriptions. This idea of energy-momentum pseudo- 
tensors was severely criticized by some people. Even it was quoted in a 
famous book [13] that: Anyone who seeks for a general formula for "local 
gravitational energy-momentum" is asking for the right answer to the wrong 
question. Misner et al. [13] showed that energy can only be localized in 
spherical systems. But later on, Cooperstock and Sarracino [H] proved that 
if energy is localizable for spherical systems, then it can be localized in any 
system. Bondi [15] argued that a non-localizable form of energy is not allowed 
in General Relativity. 

The idea of quasi- local energy was proposed by some authors [TB]-[T^. 
In this method , we can use any coordinate system while finding the quasi- 
local masses to obtain the energy-momentum of a curved spacetime. Chang 
et al. [20] proved that every energy-momentum complex can be associated 
with a particular Hamiltonian boundary term. Thus the energy-momentum 
complexes may also be considered as quasi-local. Virbhadra and his collab- 
orators [21j verified for asymptotically flat spacetimes that different energy- 
momentum complexes could give the same result for a given spacetime. They 
also found encouraging results for the case of asymptotically non-fiat space- 
times by using different energy-momentum complexes. Senovilla ^2] con- 
structed super-energy tensors for arbitrary fields in any dimension. These 
tensors had good mathematical and physical properties, and in general, the 
completely timelike component of these super-energy tensors had the math- 
ematical features of an energy density. 

It might be interesting if this problem could be explored in the alternative 
theories of gravity. Recently, some work about energy-momentum has been 
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investigated in teleparallel theory of gravity with the hope that this 

problem may be settled down in this theory. For this purpose, the teleparallel 
versions of M0ller, Bergmann, Einstein and Landau-Lifshitz prescriptions 
are derived. Sharif and Jamil [26] used these prescriptions to explore the 
energy-momentum distribution for particular spacetimes. It is concluded 
that results are consistent in some cases but no general conclusion can be 
deduced. 

The f{R) theory of gravity is another alternative theory of gravity which 
has received much attention in recent years due to its cosmologically impor- 
tant f{R) models. These models include higher order curvature invariants 
as function of Ricci scalar. It has been shown that some f{R) 

models pass solar system test. In particular, Nojiri and Odintsov [27] pro- 
posed f{R) models with negative and positive powers of the curvature. It 
is shown that the terms with positive powers of the curvature provide the 
inflationary epoch while the terms with negative powers of the curvature 
serve as effective dark energy, supporting current cosmic acceleration. They 
also discussed the consistency of some f{R) models which include the terms 
involving logarithm of scalar curvature. Cognola et al. [301 introduced a 
class of exponential f{R) models. They proved that these models passed 
all local tests, including stability of spherical body solution, non- violation of 
Newton's law, and generation of a very heavy positive mass for the additional 
scalar degree of freedom. Amendola et al. [21] derived the conditions under 
which dark energy f{R) models are cosmologically viable. Thus f{R) theory 
of gravity seems attractive due to cosmologically important f{R) models. It 
is hoped that the issue of energy-momentum localization can be addressed 
in this theory. 

In a recent paper, Multamaki et al. [22] studied energy-momentum com- 
plexes in this theory. They generalized the Landau-Lifshitz prescription 
to calculate energy-momentum in the framework of metric f{R) gravity. 
As an important special case, they evaluated the energy density for the 
Schwarzschild de Sitter spacetime. Bertolami and Sequeira [33] discussed 
some f{R) models and studied them from the point of view of the energy 
conditions and of their stability under the Dolgov-Kawasaki criterion. 

In this paper, we investigate energy distribution of some static plane sym- 
metric solutions [33] using the generalized Landau-Lifshitz energy momentum 
complex. We also explore energy density of cosmic string spacetime. These 
results are also found for some important f{R) models. The stability and con- 
stant scalar curvature conditions of these models are also discussed. The pa- 
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per is organized as follows: In section 2, we give a brief introduction about the 
field equations and the generalized Landau-Lifshitz energy-momentum com- 
plex in the context of metric f{R) gravity. In sections 3 and 4, the energy 
distribution of plane symmetric solutions and cosmic string spacetime are 
found respectively using the generalized Landau-Lifshitz energy-momentum 
complex. In the last section, we summarize and conclude the results. 

2 Generalized Landau-Lifshitz Energy- Momentum 
Complex 

The f{R) theory of gravity modifies or generalizes the general theory of 
relativity. The action for f{R) gravity is 

Here f{R) is a general function of the Ricci scalar. We note that this ac- 
tion is obtained by replacing R with f{R) in the standard Einstein-Hilbert 
action. The corresponding field equations are found by varying this action 
with respect to the metric tensor ^f^,^ 

F{R)R^,, - ^fiR)9^. - V^V,F(i?) + g^.nF{R) = kT^,, (2) 

where 

F{R) = df{R)/dR, □ = V^V^ (3) 

and represent the covariant derivative. Contracting the field equations, 
we get 

F{R)R-2f{R) + 3DF{R) = kT (4) 
and in vacuum, this reduces to 

F{R)R-2f{R) +3nF{R) = 0. (5) 

Equation ([5]) gives a relationship between f{R) and F{R). This equation 
shows that any metric with constant scalar curvature, say R = Ro, is a 
solution of the contracted equation ([5]) as long as the following equation 
holds 

F(i?o)i?o - 2/(i?o) = 0. (6) 
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This gives the condition of constant scalar curvature. For non-vacuum case, 
this is given by 

F{Ro)Ro - 2/(i?o) = kT. (7) 

These conditions are very important because these are used to check the 
acceptabihty of {R) models. This assumption of constant scalar curvature 
was firstly used by Cognola et al. [35] to investigate the solutions in f{R) 
gravity. 

The generalized Landau-Lifshitz energy-momentum complex (EMC) is 
give by [32] 

r^^ = f'iRoKl + ^(/'(^o)i?o - f{Ro))^{g'"'x' - g^'xn, (8) 

where rj^^ is the Landau-Lifshitz EMC evaluated in the framework of General 
Relativity and k = 8itG. Its 00-component turns out to be 

r°° = fiRoVA + ^{f'{Ro)Ro - /(i?o))(^^7°V + 3g''), (9) 

OK ox'- 

where is 

rTL = {-9)iT'' + tZ) (10) 
and t)e obtained from the following expression 

~ 2k '^'^ ^ aS^ Pi ^ a'l^ P5)\9 9 9 9 ) 

^ 9 9 l-L a^J- /37 + J- I3j^ aS ^ -yS^ a/3 ajJ^ -yS) 

I r,'^'^ nf^i (v^^ _L \ 

"T 9 9 aS^ IS-y ~^ ^ aS ^ -yS^ afi ^ a/3-'- -i5) 

+ rt^'(r^:,r;;,-r::,r!;,)]. (ii) 

It is mentioned here that Eq. ([8]) is the generalized formula of Landau-Lifshitz 
energy-momentum complex valid for constant scalar curvature. It would be 
worthwhile to mention here that we need cartesian coordinates to use this 
formula as some energy momentum pseudo-tensors are calculated in cartesian 
coordinates only. 

3 Energy Distribution of Plane Symmetric 
Solutions 

This section is used to evaluate energy density of some plane symmetric 
solutions found in f{R) gravity pi]. For this purpose, we use the generalized 
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Landau-Lifshitz energy-momentum complex valid for spacetimes that have 
constant scalar curvature. 



3.1 Energy Density of the 1st Solution 

The first vacuum solution (Taub's metric) is given by 

ds'^ = kix'^dt^ - dx^ - k2X^{dy'^ + dz^), (12) 

where ki and k2 are arbitrary constants. The corresponding 00-component 
takes the form 

r°° = /'(i?o)r°° + Y^(/'(^o)i?o - /(i?o))(x)i (13) 
whereas becomes 

r^L = --(^)- (14) 
Thus the 00-component of the generalized Landau-Lifshitz EMC becomes 

^ Z^f^R^) + ^(/'(i?o)i?o - /(i?o))(x)i (15) 

Now we use f{R) model to evaluate this component. It is mentioned here 
that we have some restrictions for the choice of f{R) model when R = 0. 
For example, we cannot use a model including a logarithmic function of the 
Ricci scalar and also a model which is a linear superposition of R~"^, where 
m is any positive integer. Thus we take the f{R) model [35| as follows 

f{R)=R + eR'^, (16) 

where e is any positive real number. Consequently, the 00-component of the 
generalized Landau-Lifshitz EMC reduces to 

00 ^ -5k2X^ 

3k ' ^ ^ 

Further, the stabihty condition [37], ^(i+Lfio) ~ e ^ ^' solution is 

satisfied. 
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3.2 Energy Density of the 2nd Solution 

The second vacuum solution is 

ds'' = {bx + hcfde - dx^ - e^^dy^ + dz"). 
The corresponding 00-component is 

while r^5, becomes 

Thus the 00-component of the generalized Landau-Lifshitz EMC turns out 
to be 

1 —p^'^ 1 2t 

K (x + c)^ {x-\-cy ¥[x + cY 

(20) 

For a particulat f{R) model, 

f(R) ^ R + eR'^, (21) 

this reduces to 

00 ^ ^ Y /22) 

3.3 Energy Density of the 3rd Solution 

The third solution {R ^ 0) corresponds to anti deSitter metric and is given 

by 

ds" = e2(^i-+^2)((it2 - dy' - dz") - dx\ 
Here the 00-component becomes 

= fiRoVA + ^(/'(^o)i?o - /(i?o))(|=^) (23) 



and also it follows that 



=i(-5ciV(^i^+'=^)). (24) 

K 



Thus we obtain 

= V5c,^e^(— ))/'(i?o) + limo)Ro - fm){^^)- (25) 
Now we discuss an important f{R) model given by [27] 

f{R) = R-^-hR\ (26) 
For R = Rq = 12ci^, we have 

f{Ro) = \2c^ - - 1446c/. (27) 

Inserting this value and its derivative in Eq. fl25]) . we get 

(28) 

It is mentioned here that this /(-R) model satisfy the constant scalar 
condition, i.e. F[Rq)Rq — 2f{RQ) = which implies that a = 48ci^. The 
stability condition [32j, i.e. f"{Ro) < yields 

a + b{RoY > 0. 

Since Rq = Vlc\ and a = 48ci^, it follows that 

1 + 366ci2 > 0. (29) 

Thus the model is acceptable. 

4 Energy Distribution of Cosmic String Space- 
time 

The idea of big bang suggests that universe has expanded from a hot and 
dense initial condition at some finite time in the past. It is a general cos- 
mological assumption that the universe has gone through a number of phase 
transitions at early stages of its evolution. The cosmic string spacetime has 
received serious attention in recent years due to their cosmological implica- 
tions. Here we discuss energy distribution for this cosmological model [38] in 
f{R) gravity. 
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Consider the non-static cosmic string spacetime 

ds^ = de - e^^'idp" + (1 - AGMfp'dcp' + dz'']. (30) 

We write this metric in Cartesian coordinates as it is required for the gener- 
ahzed Landau-Lifshitz EMC. In Cartesian coordinates, this becomes 

ds^ = dt' - e^'^^^^dx' - e'-'y^^^dy' 

- 1 

- e"*^^^ + 2e°*^^ ^xydxdy, (31) 

_|_ yl 

where a = 1 — 4GM with G as the gravitational constant and M as mass 

A 



per unit length of the string in the z direction and a = 2^ ^ with A as the 
cosmological constant. Also, the energy-momentum tensor is defined as 

T; = M5(x)5(y)dza(7(l, 0,0,1). (32) 

The Ricci scalar for this spacetime becomes 

Since the Ricci scalar is constant, we can find energy density of this model 
by using the generalized Landau-Lifshitz EMC. Its 00-component becomes 

= F{RMl + Uf{Ro)Ro - /(i?o))(/°). (34) 
We can evaluate by using 

rll = -9itT, + T''). (35) 
After some manipulations, it follows that 

,00 ^ 3 g 2 (^^ + qV)(^^ + Q^3:^) l^a^a l^^ " l)ix^ + a'^y^) 

(a;2 + 2/2)2 4^^^ a;2 + |/2 

x3y5(a2 _ 1)3 2xV(«' - 1)' xV(«' - l)'(a;' + aV) 



+ 



a2e"*(x2 + ?/2)5 Q2gat(^3,2 _^ ^2)5 ^2gai(2,2 _^ ^2)5 

1 ay^x^y^ix"^ + a'^y'^){y'^ + a'^x'^){a'^ - 1)^ 
2^a^ (x2 + 7/2)4 

1 a ^{x^ + a^y'^fiy'^ + a^x^)^ 
2W (a;2 + y2)4 • 



(36) 



Also, T*"^ is given by 

T^'^ = M5{x)6{y). (37) 
Using these values in Eq. (l35|) . it follows that 

x3?/5(a2 _ 1)3 2xV(«' - 1)^ ^ xV(a' - + aV) 



a2e"*(a;2 + 1/2)5 Q2gai(3,2 _^ ^2)5 gj2gat(^2,2 _^ ^2)5 

1 a 23^^2/^(2;^ + a'^y^){y^ + a2x2)(a2 - 1)2 
^ 2^a^ (x2 + 2/2)4 

- '''' 

Inserting all these values in Eg. (1511) . we have 

00 1 rprP ^.2,3atr 3 Q ^ (^^^ + « V) (Z/' + « V) 

^ = 7^-r(ii:ojae {--[-) r^— — 

2k, 4 a [x^ + y^Y 

1 ,a,2(a2-l)(x2 + a22/2) x^y\a''-lf 2x^y\a'-lf 



+ 



A^a' a;2 + |/2 a2e"*(a;2 + |/2)5 a2e°*(x2 + ?/2)5 

x'^y^{a^ - iy{x^ + a^y^) 1 a 2 (a^^ + aV)^(l/^ + a^a;^)^ 



a2e"*(x2 + ?/2)5 2' a' {x^ + y^^ 
1 ^a^^ x^y\x^ + a^y^){y^ + a^x^){a^ - If ^ ^^^^ 
+ 2^ ^^TW + MSix)Siy)} 

+ {F{Ro)Ro-f{Ro))]. (39) 



Now we discuss a well-known special case for the choice of f{R) model 

m 

f{R) =R- i-^r-'^ + i-ir-'bR"^, (40) 
where m and n are positive integers. For R = Ro = — 4A, we have 

f{Ro) = -4A + ^ - b{4Ar (41) 

and 

(4A)"+^ + an + 6m(4A)-+" 
^ ^^0) = (iA)^m • ^42) 
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Inserting these values in Eg. (15^ . it follows that 



X 



(x^ + y'^y 4 a x'^ + y'^ 

x^y^{a^ - 1)^ ^ a'^x'^y'^{x^ + a^y^){y^ + a?x^){a? - 1)^ 



+ 



x^y'^{a^ - lf{x^ + a^y^) a^{x^ + a^y^)'^{y^ + aV)^ 
a2e"*(x2 + ?/2)5 2a2(x2 + 1/2)4 

2xV(o' - 1)^ 



a2e"*(x2 + y2)5 



+ M5(x)(5(?/)} 



&(l-m)(4A)"'+"-a(l + n) 

(4A)- ^' ^ ^ 

This model must satisfy the constant curvature condition given by Eq.(I7j). 
Imposing this condition, we obtain 

a{n + 2) + b{m - 2) (4A)'^+" = (4A)"+^ - 2KM6{x)6{y) (4A)". (44) 
For a particular case, when m = 2 or 6 = 0, it reduces to 

(4A)«+i - 2KM5{x)5{y){AA)'' 



n + 2 



(45) 



It satisfies the constant scalar curvature condition which is necessary for the 
acceptability of the model. 

Now we discuss another important f{R) model [39] given by 

/(i?) = i?-aln(^) + (-l)"-^6i?". (46) 
k 

In this case, the 00-component of the generalized Landau-Lifshitz EMC takes 
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the form 

00 ^ 1 AkK-a + hkn{AKY ^ 3at r 3 a ^ (^^^ + «V)(y' + «'^') 
^ 2k^^ AkK ''''' ^ A^a' {x^ + y^f 

1 g 2 (q'-1)(^' + qV) xV(a^-l)=^ 2xV(a^-l)3 

+ ^(^). ^V(x- + aV)(y- + aV)(a2 - 1)2 ^ ^ 



+ 



+ aln(^) + ^ + 6(l-n)(4A)"]. (47) 
k k 

Also, the constant scalar curvature condition gives 
4A 

2aln(— ) - a + 4A = 6(4A)'^(n - 2) + 2KM6{x)6{y). (48) 

For n = 2 or 6 = 0, this reduces to 

2KM5{x)5{y) - 4A 



21n(^)-l 



(49) 



which satisfies the constant curvature condition necessary for the acceptabil- 
ity of the model given by (H6ll . 



5 Summary and Conclusion 

In this paper, the well-posed problem of energy-momentum localization has 
been discussed in the context of f{R) gravity. For this purpose, we use 
the generalized Landau-Lifshitz energy-momentum complex. We evaluate 
energy density of some static plane symmetric solutions by using this energy- 
momentum complex. The energy density of the cosmic string spacetime is 
also calculated. Further, this quantity is investigated for some important 
f{R) models. We have mainly considered two types of models, one with 
negative and positive powers of curvature and other including logarithmic 
term of curvature. These models have been found consistent with the solar 
system test and it has been shown that the model with negative and positive 
power of curvature unifies infiation and cosmic acceleration. The terms with 
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positive powers of curvature provide the inflationary stage while the terms 
with negative powers of curvature serves as an alternative for dark energy 
which is responsible for cosmic acceleration. 

The model with logarithm term is also cosmologically important as it 
suggests due the logarithmic term which may be responsible for the current 
acceleration of the universe. It has been shown that the chosen f{R) models 
satisfy the constant scalar curvature condition which is the necessary require- 
ment for the validity of these models. We have also explored the stability 
condition for these models. The results (fT5|) . (!20l) and (125|1 show that the 
energy density expressions are well-defined in these cases. These results can 
reduce to GR by taking f{R) = R in all the cases. We would like to men- 
tion here that we have calculated for the first time, to our knowledge, the 
energy density for a non-vacuum case (cosmic string spacetime) using the 
generalized Landau-Lifshitz energy-momentum complex. 

This work adds some knowledge about the longstanding and crucial prob- 
lem of the localization of energy. It gives the energy density expressions for 
different solutions with important f{R) models which may help at some stage 
to overcome the theoretical difficulties in the cosmological and astrophysical 
context. It would be interesting to find the Landau-Lifshitz EMC for non- 
constant scalar curvature. The extension of other EMCs in the context of 
f{R) gravity as well in other versions of f{R) gravity would also be worth- 
while. 
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